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ABSTRACT 

The  objective  of  this  analysis  is  to  discuss  the  effects  of  magnetic  field  and  its  tilted  angles  in  a  rectangular 
porous  enclosure  filled  with  nano  fluid.  The  temperature  at  moving  top  wall  is  hot  and  the  bottom  wall 
remains  with  cold  temperature,  whereas  the  other  sidewalls  are  adiabatic.  The  equations  governing  the 
considered  system  are  solved  by  the  finite  volume  method.  The  obtained  results  are  discussed  graphically  by 
using  isotherms,  streamlines  and  average  Nusselt  number.  It  shows  that  the  heat  transfer  rate  is  improved  by 
the  inclination  of  applied  magnetic  field.  In  the  shallow  cavity,  the  impact  of  Hartmann  number  is 
significant.  The  average  Nusselt  number  is  nearly  identical  when  increasing  Hartmann  number  in  the 
slender  cavity.  Regardless  of  the  permeability  and  the  aspect  ratio,  enrichment  in  the  fraction  of  solid  and 
volume  of  nano  fluid  raises  the  overall  heat  transfer  rate  in  the  occurrence  of  magnetic  field. 
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1.  INTRODUCTION 

Fluid  flow  and  heat  transfer  analysis  of  nanofluid  in  a  porous  enclosure,  is  one  of  the  most  extensively  discussed 
problems  for  the  last  few  years,  owing  to  its  uses  in  numerous  engineering  applications.  Several  researchers 
considered  the  shape  of  the  cavity  is  square  only.  Many  numerical  simulations  on  convection  of  nanofluid  in  porous 
square  cavity  have  been  discussed  by  Bourantas  etal.  [1],  Mittal  et  al.  [2],  Balia  chandra  shekar  and  Naikoti  kishan 
[3],  Kefayati  [4],  Hashemi  et  al.  [5],  Ghalambz  et  al.  [6]  and  Astanina  et  al.[7].  An  analysis  about  the  properties  of 
electrically  conducting  fluid  influenced  by  the  electromagnetic  field  is  called  Magneto-hydrodynamics  (MHD). 
Several  engineering  and  physical  applications  are  interconnected  with  the  MHD  convection.  The  orientations  of 
applied  magnetic  field  play  a  key  role  in  convectional  heat  transfer  analysis.  The  impact  of  magnetic  field  in  an 
enclosure  filled  with  fluid  or  nanofluid  or  fluid  saturated  porous  medium  has  been  analyzed  by  Chamka  [8], 
Maheswaran  etal.  [9],  Muthtamilselvan  and  Doh  [10],  Elshehabey  and  Ahmed  [11],  Hussain  [12],  Wang  etal.  [13] 
and  Grosan  etal.  [14].  Sometimes,  the  only  square  cavity  is  inadequate  for  such  heating  process,  so  the  different 
shapes  of  cavities  are  required.  For  that,  two  different  types  of  rectangular  cavity  are  considered  in  this  problem. 
This  discussion  examine  the  effects  of  Hartmann  number,  solid  volume  fraction  and  the  directions  of  applied 
magnetic  field  in  rectangular  nanofluid  saturated  porous  cavity.  Following  the  above  assessed  literature  review  and 
to  the  author’s  acquaintance,  so  far  no  such  discussion  was  carried  out. 
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2.  MATHEMATICAL  MODELLING 

Figure  1  represents  the  physical  structure  of  the  present  analysis  with  the  coordinate  system.  The  temperature  at  the  top 
wall  is  6h  while  at  the  bottom  wall  is  #csuch  that  6h  >  #f.  Adiabatic  walls  are  assumed  at  the  remaining  sides.  The  top 
wall  is  set  to  move  in  the  -y  direction  with  uniform  velocity  -£/  It  is  assumed  that  the  flow  is  steady,  Newtonian, 

incompressible  and  laminar.  The  enclosure  is  filled  with  water-Cu  nanofluid.  According  to  Boussinesq  approximation,  the 
thermal  properties  of  the  nanofluid  are  steady  excluding  the  density,  which  are  given  in  Table  1.  A  uniform  magnetic  field 

B0  is  imposed  perpendicular  to  the  adiabatic  sidewalls  and  its  tilting  angle  is  Further,  the  viscous  dissipation,  induced 

magnetic  field  effect  and  internal  heat  generation  are  taken  to  be  insignificant.  The  governing  equations  with  boundary 
conditions  for  the  considered  system  are  as  follows  [2],  [15]: 


V,  VJ 

A 


■ua 


H 


oh 

jPOTOUS  Medium^ 

Nanofluid  I 


Bn 


r 


ft 


X ,  u 

-► 


L  | 

Figure  1:  Physical  Diagram  of  the  Considered  Study. 


Table  1:  Properties  of  Nanofluid  [1],  [5] 
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u  =  v  =  0,  6  =  0C ,  y  =  0, 
u--U0,  V  =  0,  9-9h,  y  =  H, 

u  =v  =  0,  ’  x  =  0,L 

Where,  K  is  the  permeability  of  the  porous  medium,  g  is  the  gravitational  acceleration,  p  is  the  pressure,  £  is 
the  porosity  and  C p  is  the  specific  heat.  The  K  and  Forchheimer’ s  coefficient  f  are  as  follows 
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Where  dp  is  the  average  particle  size  of  the  bed.  The  non-dimensional  parameters  are  applied  in  the  equations 


(l)-(4)  are: 


X,Y=^~,  U.V=^.  T  =  4=Z-.  p- 


H 


U  n 


Oh-0,- 


^T.  H.-4 

pul  if 


/ 


/ 


V 


/ 


Pr  =  /«  =  -  2 

a.  Re2 


Thus,  (l)-(4)  becomes: 


dU  dV  „ 
dx  dY 


(5) 


dX  d  Y 


Pf  dP  +  l  Pnf 


(  32 


Pnf  dX  £  pnr  V,  Re 


2j  T  \ 


a2[/  a2t/ 
ax2  +  ar2 


/y  t/ 

PnfVf  Da  Re 


1.75  t/  (t/2+y2f  Pf  <ynS  Ha2  _7  .  ..  .  2  . 

- A  -H — - - - (V sm^cos^-t/  sm  7) 


V150  VDa 


TV  Of  Re 


(6) 


lfA+y3y 


pf  dp  +  1  /y 


/y  y  1.75  y  (i/2  +  y2)2 


3X  dY  J  pnf  dY  e  p  .  v ,  Re  pnfv ,  Da  Re  Vl50  4 Da 


'nf 

(PP)nf 


nf  v  f 
2 


+  '  •  Ri  T  +  Pf  a*  —  ( U  sin  y  cosy -V  cos 2  y) 


Pnf  P. 


f 


Pnf  Of  Re 


(7) 


ud- UV*P 


)_anf  1 

(  d 2T  d2P) 
1  _ 

j  a  f  Re  Pr 

[dx2  dY2 ) 

(8) 


with  the  boundary  conditions: 
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u  =V  =0,^L=o  X  =  0, Ar. 

dx 

U  =V  =  0,T  =  0,Y  =  0, 

U  =  - 1,  V  =  0,  T=\,  Y  =  Ar, 

Where  Re,  Gr,  Ri,  Ha ,  Pr  and  Da  stand  for  Reynolds,  Grashof,  Richardson,  Hartmann,  Prandtl  and  Darcy 
numbers,  respectively.  The  properties  of  the  nanofluid  are  [4],  [5] 

effective  density,  pnf  =  (\- %)  Pf  +  X  Ps  >  (9) 

thermal  diffusivity,  n  ^  eg  (10) 
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Based  on  Modified  Maxwell  [15]  model,  ke^  and  keq  for  spherical  nano  particles  are  expressed  by 
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where  CO  =  hnl / rs.  In  this  analysis,  the  radius  of  nano  particles (r^),  nano-layer  thickness (hnl ),  and  knl  are 
3 am,  2am  and  lOO/c^ .  The  local  (M/)  and  the  average  (Nuavg)  Nusselt  numbers  are  calculated  by  the  following  formula 
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3.  NUMERICAL  TECHNIQUE  AND  VALIDATION 


The  equations  (5)-(8)  are  solved  by  finite  volume  method  (FVM).  Coupled  pressure  and  velocity  terms  are  solved  by  Patankar’s 
[16]  SIMPLE  algorithm.  Discretised  equations  are  solved  by  TDMA  line-by-line  method.  A  grid  test  is  carried  out  for  the 

different  grid  sizes  at  Ri  -  1,  Da  =  10-2,  Ha- 50,  Pr  =  6.2,  Ar-  0.5,  y =0° ,  £*  =  0.4  and  %  =  0.06.  Table.  2  shows 

the  result  of  grid  independence  test  in  terms  of  Nuavg  and  confirms  that  the  size  of  41x81  grid  is  adequate  for  the  numerical 

calculations.  The  precision  of  the  current  numerical  calculations  evaluated  with  the  results  Arefmanesh  and  Mahmoodi  [17].  The 
results  are  shown  in  Figure  2,  which  confirms  the  accurateness  of  current  numerical  simulations. 
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Table  2:  Average  Nusselt  Number  for  Grid  Independence  Test  at  Da  =  10“2,  Ar  =  0.5,  Ha= 50, 
y=  0°,  Pr  =  6.2,  Ri  =  l,  £  =  0.4  and  ^  =  0.06. 


Grid  Size 

11  X21 

21  X41 

41  X81 

81  X  161 

0.655356 

0.603069 

0.602169 

0.602275 

Figure  2:  Comparison  of  the  Current  Results 
with  the  Results  of  Khanafer  and  Chamka  [17]. 

4.  RESULTS  AND  DISCUSSIONS 

This  article  aims  to  discussing  numerically  the  impacts  of  magnetic  field  orientations  on  combined  convection  of  nanofluid 
in  a  rectangular  porous  cavity.  The  top  wall  is  heated  uniformly  and  it  is  set  to  move  from  left  to  right.  Constant  cold 
temperature  wall  is  fixed  at  bottom.  The  vertical  sidewalls  remain  adiabatic.  The  cavity  is  filled  with  nanofluid  with 

Pr  =  6.2.  The  porosity  is  fixed  as  £*  =  0.4.  The  Ri  =  Gr/ Re 2  is  chosen  as  lthroughout  the  study.  The  results  are 
illustrated  for  various  values  of  Ha  ( 0<Ha<10 ),  solid  volume  fractions  (0  <  ;^<  0.06),  Da  (0.000 1<  Da <  0.1)  and  the 

inclination  angles  of  magnetic  field  (0°  <  y<  90°).  Moreover,  the  aspect  ratio,  Av  is  considered  as  the  ratio  of  length 
(L)  and  height  (. h ). 

Figure  3  demonstrates  the  isotherms  and  streamlines  for  various  values  of  y  at  Ha  =  50  and  Ar  -  0.5.  When 

^=0°,the  flow  inside  the  enclosure  is  characterized  by  a  counter  clockwise  (CCW)  circulating  cell  at  the  top  and  a 
clockwise  (CW)  circulating  cell  with  two  kernels  below  the  top  vortex.  The  temperature  contours  indicate  the  heat  is 
transferred  by  conduction  mode.  For  ^=30°,  the  two  kernels  of  the  bottom  cell  is  merged  and  both  the  vortices  are 

elongated  diagonally.  Further  tilting  the  magnetic  field  to  ^=60°,  the  top  vortex  increases  in  its  sizes  and  squeezes  the 
bottom  vortex.  The  isotherms  near  the  left  sidewall  at  top  part  of  the  enclosure  are  moved  towards  the  cold  wall.  When 
y  =  90° ,  the  bottom  vortex  vanishes  and  the  top  cell  engages  the  total  enclosure.  The  slope  of  isotherms  is  increased  so 
that  the  rate  of  heat  transfer  is  improved. 

Figure  4  displays  the  thermal  pattern  and  fluid  flow  for  different  inclinations  of  magnetic  field  at  Ha=  50  and 
Ar  =  2.  At  y—  0°,  the  cavity  is  filled  with  a  large  CCW  rotating  cell  below  the  small  CW  rotating  cell.  The  isotherms  are 
almost  perpendicular  to  the  vertical  wall  at  lower  part  and  slanted  type  of  isotherms  appears  at  top  part  of  the  cavity. 
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Directions  of  applied  magnetic  field  differ  from  horizontal  to  vertical,  the  fluid  flow  field  and  temperature  distribution  are 
affected.  The  bottom  cell  vanishes  and  the  top  cell  engages  the  full  cavity.  The  distribution  of  temperature  contours  is 
collapsed.  Figure  5  displays  the  heat  transfer  rate  for  various  y.  The  Nuavg  enhances  with  an  amplification  of  y.  The  heat 

transfer  rate  for  horizontal  cavity  is  higher  than  that  of  the  vertical  cavity. 


(c)  f=60° 

Figure  3:  Isotherms  and  Streamlines  for  Different 


(b)  ^=30° 


(d)  r=90° 


Y at  Ri  =  1,  Ha= 5Q  Da  =  10  2 ,  X  =  0.06  and  Ar  =  0.5. 


(a)  y=0° 


(b)  /=  3CP 


(c)  y=60° 
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(d)  y=  9(f 

Figure  4:  Isotherms  and  Streamlines  for  Different  y  Ri= \  Da  =  10  2 ,  Ha- 5Q  X  ~  0  06  and  Ar  -  2. 


Figure  5:  Average  Nusselt  Numbers  for 
Different  Ar  and  J  with  Ri - 1,  Da  =  1 0  2 
and  ^=0.06. 


Figure  6  demonstrates  the  results  of  Ha  on  streamlines  and  isotherms  at  ^=0°and  Ar—  0.5.  For  Ha  =0,  the 
enclosure  is  filled  with  two  cells,  one  below  the  other.  The  corresponding  isotherms  are  distributed  horizontally  at  lower 
part  of  the  cavity  and  skewed  type  of  isotherms  appears  at  upper  part  of  the  cavity.  When  Ha  =  20,  the  flow  is  categorized 
by  tri-cellular  flow  structure.  The  temperature  contours  are  pushed  upward.  For  Ha  =  70,  the  four  cellular  flow 
configurations  is  formed.  The  isotherms  are  straightened  and  conduction  form  of  heat  transfer  exists.  The  isotherms  and 
streamlines  for  various  Ha  with  y=0°  and  Ar—  2  are  depicted  in  Figure  7.  When  Ha  =  0,  a  large  CW  rotating  cell 
occupies  almost  the  whole  cavity  with  a  small  eddy  appears  at  the  lower  left  part  of  the  enclosure.  The  consequent 
isotherms  represent  that  high  temperature  gradient  at  lower  left  part  of  the  enclosure.  On  increasing  Ha,  the  small  vortex 
grows  bigger  and  compresses  the  primary  vortex.  The  temperature  contours  are  pushed  upward  and  extend  the  entire 
cavity. 


(a)  Ha= 0  (b)  Ha=  20 


www.tjyrc.ors 


SCOPUS  Indexed  Journal 


editor  @tjprc.  org 


1186 


S.  Muthukumar ,  Eswari  Prem  &  S.  Sureshkumar 


(c)  Ha=  50  (d)  Ha=  70 

Figure  6:  Isotherms  and  Streamlines  for  Different  Ha  at  Ri—l,  Da  =  10~2,  ^=0°,  ^  =  0.06  and  Ar= 0.5. 


(b)  Ha=  20 


(c)  Ha=  50 


(d)  Ha=10 

Figure7:  Isotherms  and  streamlines  for  different  Ha  at  Ri  =  1.  7=0°,  /)</  =  10  ’ .  ^=0.06  and 


Ar=2. 


Figure  8  illustrates  the  Nuavg  for  various  values  of  Ha.  The  heat  transfer  rate  diminishes  with  the  growing  of 
Ha.  The  effect  of  magnetic  field  is  high  for  horizontal  cavity  than  the  vertical  cavity.  The  effect  of  K  over  X  is  plotted 
in  the  Figure  9.  In  this  case,  Ha  is  chosen  as  50.  In  the  presence  of  magnetic  field,  an  enhancement  in  the  X  marks  in  the 
raise  in  Nuavg  whether  K  and  Ar  are  considered  or  not.  When  Ar  =  2,  the  higher  value  of  K  boosts  the  heat  transfer 
rate,  than  the  lower  value  of  K . 
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Figure  8:  Average  Nusselt  Numbers  for  Different  Ha  with  Ri  =  1,  Da=l(T2  and  ;^=0.06 


Figure  9:  Average  Nusselt  Numbers  for  Different  Da  and  X  with  Ri  =  1,  Ha—  50  and  ^—0°. 

5.  CONCLUSIONS 

The  consequences  of  magnetic  field  with  its  directions  in  a  rectangular  porous  cavity  filled  with  nanofluid  are  examined, 
numerically.  The  top  wall  is  heated  uniformly  and  a  constant  cold  temperature  is  maintained  at  the  bottom  wall.  The 
remaining  walls  are  fully  insulated.  The  obtained  results  are  given  below: 

•  The  inclination  angles  of  applied  magnetic  field  augment  the  heat  transfer  rate. 

•  Amplification  in  the  Hartmann  number  and  the  orientation  of  applied  magnetic  field  change  the  fluid  flow  and 
heat  transfer  characteristics. 

•  By  growing  Ha ,  the  Nuavg  is  decreased  in  the  shallow  cavity,  where  as  it  has  no  effect  in  the  slender  cavity. 

•  The  effect  of  Darcy  number  is  observed  on  heat  transfer  rate  only  when  Ar  -  2. 

•  In  spite  of  K  and  Ar,  an  enhancement  in  the  X  °f  nanofluid  raises  the  rate  of  overall  heat  transfer  in  the  existence 
of  magnetic  field. 
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